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Abstract 

^ Two different realistic F(R) modified gravity models are considered in the framework of the 

Friedmann-Lemetre- Robertson- Walker universe. The parameters of these two models are adjusted 
to reach coherence with the most recent and accurate observations of the current universe. A study 
of the growth of matter density perturbations is done, and several parametrizations of the growth 
index are developed for both models. The ansatz for the growth index given by 7 = 70 + 7i^/(l + z) 
seems to be the best parametrization for the two models considered. Finally, the values obtained 
for 70 and 71 can be used in order to characterize these two models and to differentiate them from 
others such as the Hu-Sawicki model. 
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I. INTRODUCTION 



Over the last years, several observations such as supernovae la (SNe la) [lj, large scale 
structure (LSS) [2] with baryon acoustic oscillations (BAO) [3j, cosmic microwave back- 
ground (CMB) radiation [4]-[6] and weak lensing [7J have demonstrated that our Universe is 
suffering from an acceleration in its expansion. The explanation for this accelerated expan- 
sion has become one of the most important theoretical problems for the scientific community. 
Over the last years, several suggestions have been made in attempt to solve this problem. 
One of these suggestions introduces an exotic kind of energy, called dark energy, into the 
framework of the theory of general relativity (for recent reviews in terms of dark energy, 
see [SHIP]). There is another group of theories that try to give an explanation to this late- 
time cosmic acceleration, which is based on the modification of Einstein's gravity (for recent 
reviews on modified gravity, see jTTl - fT5] ). F(R) modified gravity is among these last kinds 
of theories. 

The natural unification of inflation and late-time cosmic acceleration achieved by some 
F(R) modified gravity theories make them very attractive (see [IS]). It has also been 
demonstrated that some of these theories can be suitable candidates to explain the Universe 
we live in. In this sense, some of these F(R) gravity theories can pass the Solar System tests 
(see [T71 - fT9] ). and they can reproduce the A-cold-dark-matter (ACDM) universe in the large 
curvature regime. These theories can also accomplish the matter stability condition (see 
[T6l [201 [21]) . and they can have a stable late-time de Sitter point (see [22, 23j). Different 
observational manifestations of F(R) gravity were studied in [24J. 

As a consequence of the large number of different gravitational theories, a problem of 
distinction among some of them has appeared. The fact that different models can achieve 
the same expansion history has revealed that another tool, which may provide a way for 
discriminating among different gravitational theories, may be required. The study of the 
growth of matter density perturbations may become the tool that we need, due to the fact 
that theories with the same expansion history can have a different cosmic growth history. 
In order to characterize the growth of matter density perturbations, the so-called growth 
index 7 (see [25] ) can be very useful. 

In this work, the study of the growth history has been done for two different F(R) 
modified gravity models, and the growth index has been determined for both models. The 
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paper is organized as follows. In Sec. II, the formulations and the dynamics of modified F(R) 
gravity in the framework of a Friedmann-Lemetre-Robertson- Walker (FLRW) universe is 
briefly reviewed [26]. In Sec. Ill, two different modified F(R) gravity models are considered, 
and the values of the parameters are adjusted for both models to reach coherence with 
recent observations of the Universe. In Sec. IV, the study of the growth of matter density 
perturbations is done for these two F(R) gravity models, and several parametrizations for 
the growth index are studied for both models. Finally, a summary for this work is given in 
Sec. V. 



II. DYNAMICS OF F(R) GRAVITY IN THE FLRW UNIVERSE 

In this section, the formulation of modified F(R) gravity is reviewed in the framework of 
the FLRW universe (see [27]). The general action considered for modified F(R) gravity is 
given by 

/ = j d 4 x^ + C matt e^j , (2.1) 

M 

where Ai denotes the space-time manifold, g denotes the determinant of the metric tensor 
g^v, F(R) is a generic function of the Ricci scalar R ) £ ma tter denotes the matter Lagrangian 
and k 2 = 8ttGn = ^/M% 1 is the gravitational constant, with the Planck mass being M P \ = 
G~ 1/2 = 1.2 x 10 19 GeV. Note that, throughout this paper, the use of natural units will be 
assumed. Therefore, Kb — c — h — 1. Finally, the generic function F(R) will be given by 

F(R) = R + f(R), (2.2) 

with f(R) another generic function of R. It might be accurate here to point out that the 
R term in the expression of F(R) accounts for the Einstein-Hilbert action in the theory of 
general relativity, while the function f(R) encodes the modification of gravity. 



The variation of the action (2.1) with respect to the metric g^ v gives the field equations 
for modified F(R) gravity 



RjJLV ^ RQlLV 



1 „ 1 



2 ^ F'(R) 



l - 9lxv (F(R) - RF'(R)) + (V M V, - g^U) F'(R) 



K 2 

I rpmatter 

F'(R) ^ ' 



(2.3) 

where R^ v is the Ricci tensor, the prime denotes the derivative with respect to the Ricci scalar 
i?, V M is the covariant derivative operator associated with the metric g^ y) □(/> = g^ u V 



is the covariant d'Alembertian operator for a scalar field and T 



matter 



diag(p m , P m , P m - P m ) 



is the matter stress-energy tensor, with p m and P m being the energy density and pressure 
of matter, respectively. 



The trace of Eq.(2.3) yields 



DF'(R) 



dV P 



eff 



dF f {RY 



where 



dV e 



eff 



1 



(2F{R) - RF'(R) + K 2 T matter ) , 



(2.4) 



(2.5) 



dF'(R) 3 

with T matter = g^TJJl atter being the trace of the matter stress-energy tensor and F'(R) being 
the so-called "scalaron", i.e., effective scalar degree of freedom. 

From now on, flat FLRW space-time is assumed, with the metric given by 



ds 2 = -N(t) 2 dt 2 + a(t) 2 dx. 2 , 



(2.6) 



where N(t) is a function of the cosmic time i, which will be taken as N(t) = 1 in the 
following, and a(t) is the so-called scale factor. The Ricci scalar in this metric can be 
written as 

R = 12H 2 + 6if , (2.7) 

where H is the Hubble parameter, which is related to the scale factor a(t) through the 
expression H = a(t)/a(t), and the dot denotes the derivative with respect to the cosmic 
time t. 



From (2.3) two gravitational field equations can be obtained by considering the component 



(/x, v) = (0, 0) and the trace of (/x, v) = (with z, j = 1, 2, 3). These equations are given 
by 

3H 2 1 f 1 



K 



2H + 3H 2 



K 



F'{R) \ Pm ' 2k 2 
1 



/•"'( R) \ Pm ' 2k 2 



RF'(R) - F(R) - 6HF'(R) 
RF'(R) + F(R) + 4HF'(R) + 2F\R) 



(2.8) 
(2.9) 



By introducing the effective energy density, p e ff, and the effective pressure, P e ff, as 

f f f 



Peff 



P eff = 



1 



F'(R) 



Vr, 



F'(R) \ Pm ' 2k 2 
f 



2k 2 



RF'(R) - F(R) - 6HF'(R) 
RF'(R) + F(R) + 4HF'(R) + 2F'(R) 



(2.10) 
(2.11) 
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Eq.(2.8) can be written as 



3H 2 



2H + 3H 2 



Peff, 



effi 



(2.12) 

(2.13) 



which are the same equations that can be obtained in the framework of the theory of general 
relativity by just changing the energy density and the pressure of matter by these new 
effective energy density and effective pressure. 

In order to study the dynamics of modified F(R) gravity models in the framework of a 
flat FLRW universe, it may be useful to introduce the new function yn(z) given by 



Pm(o) m 2 



(2.14) 



with z = l/a(t) — 1 being the redshift, p m (o) being the matter energy density at the present 
time, m 2 being the mass scale given by 

~2, 



~ 2 K Pm(0) 

m = — 



~ 1.5 x 10~ 67 eV 2 



and, finally, x being defined as 



X 



_ Pr(0) 



Pm(0) 



3.1 X 10" 



(2.15) 



(2.16) 



where p r ( ) is the radiation energy density at the present cosmic time. 



Introducing Eq.(2.14) into the first equation in (2.12) yields 

d 2 y H (z) , j dy H (z) 



dz 2 



+ Ji 



where 



1 



+ 
1 



dz 



1 



+ -h (vh(z)) + J 3 = 0, 



1 - F'(R) 



y H + (z + l) 3 + x(z + l) 4 Qm 2 F"(R)_ 
1 2 - F'(R) 



Ji 

32 (z + l) 2 [y H + {z + If + X (z + l) 4 3m 2 F"(R)_ 
J 3 = -3(z + 1) 

(1 - F'(R))((z + l) 3 + 2 X (z + l) 4 ) + (R- F{R))/{3fh 2 ) 

(^ + i) 2 (^ + (^ + i) 3 + x(^ + i) 4 ) 

Furthermore, the Ricci scalar is expressed as 



Qfn?F"{R) ' 



R = 3m 



^2 



4y H (z)-(z + l) d ^l + (z + lf 



(2.17) 

(2.18) 
(2.19) 

(2.20) 
(2.21) 
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In deriving this equation, we have used the fact that — [z + l)H(z)d/dz = 
H(t)d/d(\na(t)) = d/dt, where H could be an explicit function of the redshift as H = H(z), 
or an explicit function of the time as H — H(t). 



Once the function yn(z) is obtained from Eq.(2.17), it can be used to calculate the 
equation of state of dark energy uu DE through the expression given by 

uj de (z) = = -1 + -(z + 1 — — , 2.22 

Pde 3 y H (z) d(z) 

and the dark energy density parameter Q,be given by 

toDE{z) = ^ = — f (2.23) 
which are two essential quantities in determining whether a theory is realistic or not. 



III. REALISTIC F(R) MODELS 

In this section, two different kinds of modified F(R) gravity models will be considered. 
The parameters of these models will be set in order to reproduce recent observations of our 
current Universe. 

In j6], Komatsu et al determined important cosmological parameters by combining 
the seven-year WMAP data with the latest distance measurements from the (BAO) in the 
distribution of galaxies, the Hubble constant (H ) measurement and the last observations 
coming from the luminosity distances out to high-z type la supernovae (SN). The determined 
values for the dark energy equation of state parameter uj de and for the dark energy density 
parameter Vt DE are given by 

u DE = -0.980 ± 0.053 from (WMAP+BAO+SN) , 
n DE = 0.725 ± 0.016 from (WMAP+BAO+H ) . (3.1) 
From now on, these results will be used as a constraint for the two model parameters. 



A. First F(R) model 

In the first place, a model that appeared in [28] which could unify inflation and current 
acceleration will be considered. This model is given by 
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where a, /3 and 7 are positive constants and m, n and / are positive integers satisfying the 



condition m + I > n. The model given by (3.2) is a generalization of the Hu-Sawicki model 
(see [29J) which has been proposed by Hu and Sawicki as a model which is in agreement 



with the constraints imposed by the Solar System tests. Model (3.2) can be reparametrized 
by choosing n = /, /3 = 2A/(6A) n and 7 = l/(bk) n yielding 



F(R) = R - 2A 1 - 



1 



+ 



aR 



m+n 



(3.3) 



i + (fi)V i + (fi)"' 

with A being the current cosmological constant. It is worth noting that the new constant b 
can be negative when n is a positive even; in any other it must be positive. 



The next step should be to solve Eq.(2.17) for model (3.3) and to find out the constraints 



on the values of the model parameters needed to fulfill the conditions given by (JTTTj) . Un- 



fortunately, because Eq.(2.17) cannot be solved in an analytical way for the F(R) model 



given by (3.3), this is not possible. Thus, the way to solve this problem is to suggest a 



set of parameters for (3.3), to solve Eq.(2.17) numerically and to check if the results are in 



accordance with (III). 



For the model given by (3.3), the following values for the parameters have been chosen: 

3 



n = 4, m — 1, b = 



a = 10 



- 10 m- 8 , 



(3.4) 



with A = 7.93m 2 in accordance with [6]. In order to obtain the initial conditions needed to 



solve Eq.(2.17) numerically, we may evaluate the dark energy density Pde = Peff — Pm from 
Eq.(2.10) at the matter dominated era (high redshifts) by putting R = 3m 2 (l + z) 3 . In the 



case of the first model with the set of parameters given by (3.4), the initial conditions can 
be written as 



Vh\z) 



dy H (z) 



A 

3m 2 



-81mfi 4 (l + z) c 



dz 



= -243mfi 4 (l + z) 2 . 



(3.5) 



In the case of model (3.3) with the set of parameters given by (3.4) and the initial conditions 



given by (III A), I set Zi = 3.40, obtaining ^^(0) = —1.000 and Qde — 0.725, which are in 



accordance with the observational data given by flllip. Note that it is hard to solve Eq.(2.17) 



for higher values of the redshift due to the large frequency of the dark energy oscillations. 



In the following, model (2.2) with the set of parameters given by (3.7) will be called 
model I. 
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B. Second F(R) model 



The second model considered |30j is given by 
F(R) = R-a 



-a j 



tanh 



bi(R-Ri) 



tanh 



biRi 



(3.6) 



It will be assumed that Ri S> Rq, aj 3> cto and 6/ <C 6q, with S> 1. By choosing 



2A = a [1 + tanh (^f 41 )] and 2A 7 = a/ [l + tanh (^f 1 )] , Eq.Q reduces to 

1 



F(R) = R-2A 



1 



_ tanh (^zM 
1 + tanh pfaj" 



2A 7 



1 



1 _ tanh (M^i) 
1 + tanh (^) 



(3.7) 



where A is the current cosmological constant, while A/ accounts for the effective cosmological 
constant in the early Universe. In the following it will not be considered the last part, the 



one that accounts for inflation, in (3.7). Thus, the second model we take into consideration 
will be the one given by 



F(R) = R-2A 



1 



1 _ tanh (M^o) 
1 + tanh (^fa) 



(3i 



As already shown in the first model, analytical solutions for Eq.(2.17) cannot be found 



for (3.8). The procedure followed in order to solve the problem is the same one used in the 



previous subsection; i.e., a set of parameters will be chosen for model (3.8), then Eq.(2.17) 



will be solved numerically, and, finally, I will check whether the results are in accordance 



with (III) or not 



For the model (3.8), I set 



R = l(T 66 ei/ 2 , b = 1.WRq\ A = 7.93m 2 



(3.9) 



Following the same steps as in the previous subsection, it is found that the initial conditions 
are given by 

A 



1 — tanh ( b 



Vh(z)\. 
dy H (z) 



dz 



3fh 2 
3bA 



1 



1 + tanh (*fa) 



cosh ( 6 aft a (i+«)'-*b 



1 + tanh (*fa) 



-(1 + ,) 2 . 



(3.10) 
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And, finally, for model (3.8) with the set of parameters (3.9) and initial conditions given 



by (III B), I set Zi = 2.51, obtaining ^^(0) = —0.969 and CIde = 0.735, which are also in 



accordance with the observational data given by (III). 



From now on, model (2.5) with the set of parameters given by (3.9) will be called model 



II. 



IV. GROWTH OF MATTER PERTURBATIONS: GROWTH RATE 

In this section, the growth of matter density perturbations, 8 = for model I and 
model II is studied. Since it is known that many different gravitational theories can mimic 
the ACDM unverse, which is commonly accepted as the Universe in which we live, the study 
of the growth history of these theories may be considered an essential tool to discriminate 
among them. 

Under the subhorizon approximation, the matter density perturbation 5 = — satisfies 
the following equation [3]] (and references therein): 

6 + 2H6 - 47rG e ff(a, k)p m 6 = 0, (4.1) 

with k being the comoving wave number and G e ff(a, k) being the effective gravitational 
"constant" given by 

(k 2 /a 2 ) (F»(R)/F'(R)) 



G 

G e $(a,k) 



F>(R) 



1 + 



(4.2) 



l + 3(k 2 /a 2 )(F"(R)/F'(R))_ ' 
In Figs. [l]and[2j the cosmological evolution of the ratio G e ff/G as a function of redshift 
z and the comoving wave number k for both model I and model II is depicted. 

The appearance of the comoving wave number k in the expression of the effective gravi- 
tational constant G e ff has a huge importance due to the fact that now the evolution of the 
matter density perturbations also depends on k. This kind of dependence does not appear 
in the framework of general relativity. This fact can be easily checked by taking F(R) = R 



in Eq. (4.2). 



In deriving Eq. (4.1), I assume the subhorizon approximation (see [32]). for which the 
comoving wavelengths A = a/k are considered to be much shorter than the Hubble radius 
H~ l . In terms of the comoving wave number, the subhorizon approximation can be written 
as follows: 



% » H 2 . (4.3) 
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FIG. 1: Contour plot of the effective gravitational constant G e ^/G as a function of z and 
logfc(Mpc _1 ) for model I (a) and model II (b). 
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FIG. 2: Cosmological evolution of G e ^/G as a function of z for k = O.lMpc 1 (blue line), k — 
O.OlMpc" 1 (red line) and k = O.OOlMpc" 1 ( green line) for model I (a) and model II (b). 

For model I, the subhorizon approximation states that k 3> 0.000116Mpc -1 . For model 



II, it states that k 3> 0.000118Mpc 1 . In order to fulfill (4.3), the values considered for 



k in this work will always satisfy \ogk > —3, with k written in Mpc -1 , for both model I 
and model II. From now on, the expression log k (Mpc -1 ) will be used to specify taking the 
logarithm of fc, with k written in Mpc -1 . On the other hand, deviations from the linear 
regime have to be taken into account [33] when log£;(Mpc -1 ) > —1. Thus, the range of 
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(b) / g (model II) 



FIG. 3: Contour plot of the growth rate f g as a function of z and log/c(Mpc 1 ) for model I (a) 
and model II (b). 
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FIG. 4: Cosmological evolution of the growth rate f g as a function of z for k = O.lMpc 1 (blue 
line), k = O.OlMpc" 1 (red line) and k = O.OOlMpc" 1 ( green line) for model I (a) and model II (b). 

values considered for k throughout this work for both model I and model II is 

-3 < logfc(Mpc _1 ) < -1. 



Equation (4.1) can be written in a different way by using the so-called growth rate f g 



given by f g = din 5 /din a. In terms of this growth rate, Eq. (4.1) reduces to 



df g (z) + ( l + zdH(z) 2 



dz 



H(z) dz 



f g (z) 3fh 2 (l + z) 2 G eS (a(z),k) 



1 + z + 2 



H 2 (z) 



G 



0. 



(4.4) 
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FIG. 5: Constant growth index as a function of log /c(Mpc _1 ) for model I (a) and for model II (b). 
The bars express the 68% CL. 

This equation can be solved numerically for model I and model II by imposing the con- 
dition that at high redshift the results for the ACDM universe are recovered. In Figs. [3] and 
[4j the growth rate is shown as a function of the redshift z and the comoving wave number 
k for model I and model II. 

The next step should be to use the growth rate of these two models to compare them, 



but a new problem comes up when we face this task. Equation (4.4) usually must be solved 



numerically because of its complexity. This means that, generally, we will not have an 
analytic expression for the growth rate to deal with. In order to compare and discriminate 
among different theories it would be helpful to have one or more parameters that characterize 
their growth history. 



V. CHARACTERIZING THE GROWTH HISTORY: GROWTH INDEX 

In this section the concept of the growth index is developed. The growth index 7 appears 
as an important quantity in characterizing the growth of matter density perturbations. 

In order to compare the growth of matter density perturbations between different theories, 
the so-called growth index 7 appears. This index is given by 

f g (z) = Q m (z)^\ (5.1) 
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FIG. 6: Cosmological evolutions of the growth rate f g (red line) and Qja (blue line) with 7 = 70 as 
functions of the redshift z in model I for k = O.lMpc" 1 (a), k = O.OlMpc" 1 (c) and k = O.OOlMpc" 1 
(e), and those in model II for k = O.lMpc" 1 (b), k = O.OlMpc" 1 (d) and k = O.OOlMpc" 1 (f). 

where ^ m (^) = 8 g^ m is the matter density parameter. The growth index 7 cannot be 
directly observed, but it could have a huge importance in discriminating among different 
gravitational theories; it can be inferred from the observational data of both the growth 
factor f g (z) and the matter density parameter Q m (z) at the same redshift z. 

As it was done in [27] , different parametrizations for the growth index 7 will be considered 
for both model I and model II. In a first stage, a constant 7 will be assumed jSH 135] : 
afterwards, a linear dependence [36j given by j(z) = 70 + % • z will be considered; and, 
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(a) model I (b) model II 



FIG. 7: Cosmological evolution of the relative difference e r = 1 g ^ m| with 7 = 70 for k = O.lMpc 
(blue line), k = O.OlMpc" 1 (red line) and k = O.OOlMpc" 1 ( green line) in model I (a) and model 

11 (b). 

finally, an ansatz of the type j(z) = 70 + 71 • z/(l + z) will be suggested. 

In what follows, we will study these different par ametrizat ions of the growth index for 
model I and model II. 



A. 7 = 70 

We consider the ansatz for the growth index given by 

7 = 7o, (5.2) 

where 70 is a constant. 



The results obtained by fitting Eq. (5.1) to the solution of Eq. (4.4) for different values 
of the comoving wave number k for model I and model II are shown in Fig. [5j where the 
points denote the median value while the bars express the 68% confidence level (CL). It 
can be easily observed that both models exhibit a strong and quite similar dependence on 
log fc(Mpc _1 ). Moreover, 7 seems to be worse determined for model II. 

The cosmological evolutions of the growth rate f g (z) and the expression ^ m (z) 70 as func- 
tions of the redshift z for several values of the comoving wave number k for model I and 
model II are depicted in Fig. [6| From what is shown in this figure it is clear that the worst 
fit is given for the highest value of the comoving wave number k. 
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In order to clarify the results obtained, the relative difference between f g (z) and ^ m (z) 70 
is defined as 

«,(».*) - l/ ' ( '-* ) ry* )T| . (5.3) 

In Fig. [7| the cosmological evolution of e r as a function of z for the same values of k in both 
models are shown. For model I the relative difference is less than 13% for log fc(Mpc _1 ) = 
-1, 3.5% for logfc(Mpc _1 ) = -2 and 2% for logfc(Mpc _1 ) = -3; while for model II, the 
highest value of e r is 14% for log fc(Mpc _1 ) = —1, 7% for log fc(Mpc _1 ) = —2 and 3% for 
log£;(Mpc -1 ) = —3. Thus, two points can be made. First of all, the fits for model I are, 
generally, better than the ones for model II. Second of all, the fits are better for lower values 
of log fc(Mpc _1 ) for both models. 

B. 7 = 70 + 71 ' z 

In this subsection, the case of a growth index given by 

7 = 7o + 7i ' ^ ( 5 - 4 ) 

where 70 and 71 are constants, will be studied following the same steps taken in the case of 
a constant growth index. The results obtained with this ansatz should improve those for a 
constant growth index. 

In Fig. [sj the parameters 70 and 71 for several values of log/c(Mpc _1 ) in both models 
are shown. For model I, 70 exhibits a clear dependence on log /c(Mpc _1 ), while 71 is almost 
constant for —3 < log fc(Mpc _1 ) < —1.75. For model II, the dependence on log fc(Mpc _1 ) 
is strong for both 70 and 71 throughout the range of values considered for log fc(Mpc _1 ). 
It may also be noted that the parameter 71 gives the main difference between model I and 
model II. 

In Fig. [9J the cosmological evolutions of the growth rate f g (z) and ^ m (^) 7 ^ as functions 
of the redshift z together for model I and model II are depicted. Compared with the fits for 
a constant growth index, it can be easily noticed that the linear ansatz improves the results 
obtained, especially in the case of log^Mpc- 1 ) = -1. 

The cosmological evolution of the relative difference e r as a function of z for several values 



of k in model I and model II is shown in Fig. 10, For model I the relative difference is less 



than 3.5% for log/c(Mpc 1 ) = —1 and 1.5% for log/c(Mpc 1 ) < —2; while for model II, the 
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FIG. 8: Growth index fitting parameters in the case 7 = 70 + 71 • z as a function of log A:(Mpc _1 ) 
for model I [(a) and (c)] and model II [(b) and (d)]. The legend is the same as Fig. [5j 

highest value of e r is 2.5% for logfc(Mpc _1 ) = -1, 1.5% for logfc(Mpc _1 ) = -2 and 2.5% for 
log /c(Mpc _1 ) = —3. It might be accurate to note that these results improve those obtained 
for the case given by 7 = 70, particulary those for the case log/c(Mpc _1 ) = —1. 



C. 7 = 70 + 71 



l+z 



Finally, we assume the following ansatz for the growth index: 



7 = 7o + 7i 



(5.5) 
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FIG. 9: Cosmological evolutions of the growth rate f g (red line) and QZi (blue line) with 7 = 
7o + 7i • z as functions of the redshift z in model I for k = O.lMpc" 1 (a), k = O.OlMpc" 1 (c) 
and k = O.OOlMpc" 1 (e), and those in model II for k = O.lMpc" 1 (b), k = O.OlMpc" 1 (d) and 
k = O.OOlMpc- 1 (f). 

where 70 and 71 are constants. 

The parameters 70 and 71 for several values of log k are shown in Fig. [TT] for both mod- 
els. For model I, as it happened in the linear case, 70 exhibits a strong dependence on 
log /c(Mpc _1 ), while 71 is almost constant for —3 < log fc(Mpc _1 ) < —2. In the case of 
model II, the dependence on log/c(Mpc _1 ) is clear for both parameters 70 and 71. 

The cosmological evolutions of the growth rate f g (z) and ^ m (z) 7 ^ in model I and model 



17 



0.5 1.0 1.5 2.0 2.5 3.0 



0.5 1.0 1.5 2.0 
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(a) model I 



FIG. 10: Cosmological evolution of the relative difference e r 



(b) model II 

fg — 



-j with 7 = 70 + 71 • z for 

k = O.lMpc -1 (blue line), k = O.OlMpc -1 (red line) and k = O.OOlMpc -1 (green line) in model I 
(a) and model II (b). 



II for several values of k are depicted in Fig. [12} The fits in this case improve the results 
obtained for a constant growth index, but they seem similar to those obtained for the linear 
case. 

As in the previous subsections, the relative difference e r for several values of k in 



model I and model II is shown in Fig. 13 in order to analyze the fits quantitatively. 
For model I the relative difference is less than 3% for log fc(Mpc _1 ) = —1 and 1.5% for 
log /c(Mpc _1 ) < —2; while for model II, the highest value of e r is 5% for log fc(Mpc _1 ) = — 1, 
3% for log fc(Mpc _1 ) = —2 and 3.5% for log/c(Mpc _1 ) = —3. These results improve those 
obtained for the case given by 7 = 70. For model I the results are very similar to the 
corresponding ones of 7 = 70 + 71 • z, but for model II they are worse. 

To conclude, it is important to point out that three parametrizations for the growth index 
have been studied for both model I and model II. As it may have been expected, the fits 
obtained for a constant growth index give the worst results. The results obtained for the 
other two ansatz considered, i.e. 7 = 70 + 71 • z and 7 = 70 + 7i • j^, are quite similar 
for model I, but 7 = 70 + 71 • z gives better fits for model II than those corresponding to 
the ansatz 7 = 70 + 7i • In conclusion, the linear ansatz, 7 = 70 + 71 • z, is the best 
parametrization for the growth index for the two models considered. 
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FIG. 11: Growth index fitting parameters in the case 7 = 70 + 71 • as a function of log /c(Mpc _1 ) 
for model I [(a) and (c)] and model II [(b) and (d)]. The legend is the same as Fig. [5j 

VI. DISCUSSION 

In this section, the content of the paper is summarized and the results obtained for model 
I and model II are analyzed. 



Two models of F(R) modified gravity given by (3.3) and (3.8) have been considered 
throughout this work. The parameters of these models have been set in order to agree with 



current observational data coming from [6j. Model (3.3) with the set of parameters given 



by (3.4) is the so-called model I, while model (3.8) with the set of parameters given by (3.9) 



is the so-called model II. 
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FIG. 12: Cosmological evolutions of the growth rate f g (red line) and ftZi (blue line) with 7 = 
7o + 7i • 1^ as functions of the redshift z in model I for k = O.lMpc -1 (a), k = O.OlMpc -1 (c) 
and k = O.OOlMpc" 1 (e), and those in model II for k = O.lMpc" 1 (b), k = O.OlMpc" 1 (d) and 
k = O.OOlMpc- 1 (f). 

The growth of matter density perturbations has been studied for model I and model II. 
The so-called growth rate has been obtained numerically for both models and three ansatz 



for the so-called growth index have been considered. In Figs. [5j |8| and [Tit the results of the 
different parametrizations for the growth index are shown for both models. 



To determine which ansatz of those considered for the growth index fits better Eq. (5.1) 



to the solution of Eq. (4.4), the results obtained for the three parametrizations have been 
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FIG. 13: Cosmological evolution of the relative difference e r 



(b) model II 



with 7 = 70 + 7i • ^ 



for 



k = O.lMpc -1 (blue line), k = O.OlMpc -1 (red line) and k = O.OOlMpc -1 (green line) in model I 
(a) and model II (b). 

analyzed in the previous section. The ansatz given by 7 = 70 + 71 • z seems to be the best 
choice for both models. 

Thus, in order to discriminate between model I and model II (or with the models con- 
sidered in [27] ) using the growth history, the values of 70 and 71 in 7 = 70 + 71 • z could 
have an essential importance. In Fig. [14] these two parameters, 70 and 71, are depicted for 
the two models together. We see that the behavior of 70 is very similar for both models. 
However, the values for 71 are totally different for model I to model II and they could be 
used to discriminate between these two models. 



One final note must be made. As it was pointed out in Eq.(4.1), the evolution of matter 



density perturbations for F(R) modified gravity theories depends on the comoving wave 
number fc, which does not occur in the framework of general relativity. Throughout this 
work, this fact has been confirmed, in the first place with the results obtained for the growth 
rate / g , and finally, with those obtained for each of the three parametrizations considered for 
the growth index 7. Nevertheless, these parametrizations do not depend on the comoving 
wave number k ) and it may be very interesting in the future to propose some scale-dependent 
ansatz for these F{R) modified gravity theories. 



21 



0.55 



0.50 



0.45 



0.40 



Model I 
Model II 



-0.05 



-0.10 



-0.15 



-0.20 



* i * 



* * 4 



Model I 
Model II 



i 



-3.0 -2.5 -2.0 -1.5 -1.0 
Log k(Mpc _1 ) 



-3.0 -2.5 -2.0 -1.5 -1.0 

Log k(Mpc _1 ) 



(a) (b) 

FIG. 14: Growth index fitting parameters in the case 7 = 70 + 71 • z as a function of log &;(Mpc _1 ). 
The legend is the same as Fig. [5j 

Acknowledgments 

I would like to thank the referee for comments and criticisms that led to the improvement 
of a first version. I would also like to thank Sergei Odintsov for suggesting the problem and 
for giving me some ideas to carry out this task. I would like to acknowledge a JAE fellowship 
from CSIC. 



[1] S. Perlmutter et al [SNCP Collaboration], Astrophys. J. 517, 565 (1999) | arXiv:astn>l 

|ph/9812133 |; A. G. Riess et al [Supernova Search Team Collaboration], Astron. J. 116, 1009 

(1998) |arXiv:astro-ph/98052011 . 
[2] M. Tegmark et al [SDSS Collaboration], Phys. Rev. D 69, 103501 (2004) [arXiv:astro- 

ph/0310723 ; U. Seljak et al [SDSS Collaboration], Phys. Rev. D 71, 103515 (2005) 

|arXiv:astro-ph/0407372] . 
[3] D. J. Eisenstein et al [SDSS Collaboration], Astrophys. J. 633, 560 (2005) | arXiv:astro- 

|ph/050117l] . 

[4] D. N. Spergel et al [WMAP Collaboration], Astrophys. J. Suppl. 148, 175 (2003) | arXiv:astro- 

ph/0302209|; Astrophys. J. Suppl. 170, 377 (2007) | a rXiv:astro-ph/0603449l . 
[5] E. Komatsu et al [WMAP Collaboration], Astrophys. J. Suppl. 180, 330 (2009) 



22 



[arXiv:0803.0547| [astro-ph]]. 
[6] E. Komatsu et al [WMAP Collaboration], Astrophys. J. Suppl. 192 (2011) 18 

[arXiv: 1001. 4538] [astro-ph.CO]]. 
[7] B. Jain and A. Taylor, Phys. Rev. Lett. 91, 141302 (2003) |ar Xiv:astro^p h/0306046|. 
[8] M. Li, X. D. Li, S. Wang and Y. Wang, Commun. Theor. Phys. 56, 525 (2011) |arXi v:1103.5"870 

[astro-ph.CO]]. 

M. Kunz, |aIXiv:1204.5482| [astro-ph.CO]. 

K. Bamba, S. Capozziello, S. Nojiri and S. D. Odintsov, arXiv:1205.3421 [gr-qc]. 



[9 
[10 
[11 

[12; 
[13 
[14 

[is: 

[16 



[17: 
[is: 

[19 
[20 

[21 
[22 



S. Nojiri and S. D. Odintsov, Phys. Rept. 505, 59 (2011) [arXiv: 101 1.0544 [gr-qc]]; eConf 



C0602061, 06 (2006) [Int. J. Geom. Meth. Mod. Phys. 4, 115 (2007)] |arXiv:hep-th/0601213|. 
S. Capozziello and V. Faraoni, Beyond Einstein Gravity: A Survey of Gravitational Theories 
for Cosmology and Astrophysics (Springer, Berlin, 2010). 

T. Clifton, P. G. Ferreira, A. Padilla and C. Skordis, Phys. Rept. 513, 1 (2012) 



| arXiv: 1 106 . 2476 [astro-ph. CO] ] . 



S. Capozziello and M. De Laurentis, Phys. Rept. 509, 167 (2011) |arXiv:1108.6266 [gr-qc]]. 
S. Capozziello, M. De Laurentis and S. D. Odintsov, |arXiv: 12 06.4842 [gr-qc]. 
S. Nojiri and S. D. Odintsov, Phys. Lett. B 657 (2007) 238 | arXiv:0707.1941| [hep-th]] ; S. Nojiri 
and S. D. Odintsov, Phys. Rev. D 77 (2008) 026007 |arXi^0710l738] [hep-th]], G. Cognola, 
E. Elizalde, S. Nojiri, S. D. Odintsov, L. Sebastiani and S. Zerbini, Phys. Rev. D 77 (2008) 



046009 |arXi v:0712.4017 [hep-th]]. 

S. Nojiri and S. D. Odintsov, Phys. Rev. D 68 (2003) 123512 |hep- th/0307288| ; 
T. Chiba, Phys. Lett. B 575, 1 (2003) |arXiv:astro-ph/0307338|. 



T. Chiba, T. L. Smith and A. L. Erickcek, Phys. Rev. D 75, 124014 (2007) |arX iv:astro- 



ph/0611867|; G. J. Olmo, Phys. Rev. D 75, 023511 (2007) |gr-qc/0612047|. 

A. D. Dolgov and M. Kawasaki, Phys. Lett. B 573, 1 (2003) [arXiv:astro-ph/0307285l ; 



V. Faraoni, Phys. Rev. D 74, 104017 (2006) |arXiv:astro-ph/0610734|. 



Y. S. Song, W. Hu and I. Sawicki, Phys. Rev. D 75, 044004 (2007) |arXiv:astro-ph/0610532|. 
V. Muller, H. J. Schmidt and A. A. Starobinsky, Phys. Lett. B 202, 198 (1988); V. Faraoni 



and S. Nadeau, Phys. Rev. D 72, 124005 (2005) arXiv :gr-qc/051 1094|. 
[23] G. Cognola, E. Elizalde, S. D. Odintsov, P. Tretyakov and S. Zerbini, Phys. Rev. D 79 (2009) 
044001 |arXi v:0810.4989l [gr-qc]]. 



23 



[24] A. de la Cruz-Dombriz, A. Dobado, and A. L. Maroto, Phys. Rev. D 77, 123515 (2008); L. 
Pogosian, and A. Silvestri, Phys. Rev. D 77, 023503 (2008); idem, Phys. Rev. D 81, 049901(E) 
(2010); H. Oyaizu, Phys. Rev. D 78, 123523 (2008); H. Oyaizu, M. Lima, and W. Hu, Phys. 
Rev. D 78, 123524 (2008); F. Schmidt, M. Lima, H. Oyaizu, and W. Hu, Phys. Rev. D 79, 
083518 (2009); G. B. Zhao, B. Li, and K. Koyama, Phys. Rev. D 83, 044007 (2011); K. Koyama, 
A. Taruya, and T. Hiramatsu, Phys. Rev. D 79, 123512 (2009); A. Taruya, T. Nishimichi, S. 
Saito, and T. Hiramatsu, Phys. Rev. D 80, 123503 (2009); F. Schmidt, Phys. Rev. D 78, 
043002 (2008); G. B. Zhao, L. Pogosian, A. Silvestri, and J. Zylberberg, Phys. Rev. D 79, 
083513 (2009); F. Schmidt, A. Vikhlinin, and W. Hu, Phys. Rev. D 80, 083505 (2009); A. 
Borisov, and B. Jain, Phys. Rev. D 79, 103506 (2009); S. Ferraro, F. Schmidt, W. Hu, Phys. 
Rev. D 83, 063503 (2011); K. W. Masui, F. Schmidt,, U. L. Pen, and P. McDonald, Phys. Rev. 
D 81, 062001 (2011); T. Giannantonio, M. Martinelli, A. Silvestri, and A. Melchiorri, JCAP 

04, 030 (2010); K. Yamamoto, G. Nakamura, G. Hutsi, T. Narikawa, and T. Sato, Phys. Rev. 
D 81, 103517 (2010); E. Beynon, D. J. Bacon, and K. Koyama, Mon. Not. Roy. Astron. Soc. 
403, 353 (2010); A. M. Nzioki, P. K. S. Dunsby, R. Goswami, and S. Carloni, Phys. Rev. D 
83, 024030 (2011); X. Fu, P. Wu, and H. Yu, Eur. Phys. J. C 68, 271 (2010); T. Narikawa, and 
K. Yamamoto, Phys. Rev. D 81, 043528 (2010); ibidem, Phys. Rev. D 81, 129903(E) (2010); 

5. A. Thomas, S. A. Appleby, and J. Weller, JCAP 03, 036 (2011); L. Lombriser, F. Schmidt, 



T. Baldauf, R. Mandelbaum, U. Seljak, and R. E. Smith, arXi v:1111.2 020; S. Camera, A. 
Diaferio, and V. F. Cardone, JCAP 07, 016 (2011) 
[25] E. V. Linder, Phys. Rev. D 72, 043529 (2005) lastro-ph/0507263l . 

[26] S. Capozziello, S. Nojiri and S. D. Odintsov, Phys. Lett. B634, 93 (2006), |hqTth/0512118| 
S. Capozziello, S. Nojiri, S. D. Odintsov and A. Troisi, Phys. Lett. B639, 135 (2006) , |astrc> 
ph/0604431 



[27] K. Bamba, A. Lopez-Revelles, R. Myrzakulov, S. D. Odintsov and L. Sebastiani, 



arXiv: 1207. 1009 [gr-qc]. 



[28] S. Nojiri and S. D. Odintsov, Phys. Rev. D 77 (2008) 026007 |arXiv: 07iai738| [hep-th]]. 
[29] W. Hu and I. Sawicki, Phys. Rev. D 76 (2007) 064004 |arXiv:0705.1158| [astro-ph]]. 



[30] S. Nojiri and S. D. Odintsov, Prog. Theor. Phys. Suppl. 190 (2011) 155 [arXiv: 1008.4275 
[hep-th]]. 

[31] S. Tsujikawa, Phys. Rev. D 76 (2007) 023514 | arXiv:0705. 10321 [astro-ph]]. 

24 



[32] G. Esposito-Farese and D. Polarski, Phys. Rev. D 63 (2001) 063504 |arXiv:gr-qc/00 09034|. 
[33] V. F. Cardone, S. Camera and A. Diaferio, JCAP 1202 (2012) 030 |i arXiv:1201.3"272] [astro- 
ph]]. 

[34] N. Kaiser, Mon. Not. Roy. Astron. Soc. 227 (1987) 1. 
[35] A. J. S. Hamilton, astro-ph/9708102. 



[36] D. Polarski and R. Gannouji, Phys. Lett. B 660 (2008) 439 arXiv:0710.1510 [astro-ph]]. 



25 



